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We study procedures for the optimization of efficiency of light storage and retrieval based on the 
dynamic form of electromagnetically induced transparency (EIT) in warm Rb vapor. We present a 
detailed analysis of two recently demonstrated optimization protocols: a time-reversal-based itera- 
tion procedure, which finds the optimal input signal pulse shape for any given control field, and a 
procedure based on the calculation of an optimal control field for any given signal pulse shape. We 
verify that the two procedures are consistent with each other, and that they both independently 
achieve the maximum memory efficiency for any given optical depth. We observe good agreement 
with theoretical predictions for moderate optical depths (< 25), while at higher optical depths the 
experimental efficiency falls below the theoretically predicted values. We identify possible effects 
responsible for this reduction in memory efficiency. 

PACS numbers: 42.50.Gy, 32.70.Jz, 42.50.Md 



I. INTRODUCTION 

The ability to store light pulses in matter and then 
retrieve them while preserving their quantum state is an 
important step in the realization of quantum networks 
and certain quantum cryptography protocols [l|,|2|. Map- 
ping quantum states of light onto an ensemble of identical 
radiators {e.g., atoms, ions, solid-state emitters, etc.) of- 
fers a promising approach to the practical realization of 
quantum memory [S EL IS HI • Recent realizations^fstor- 
age and retrieval of single-photon wave packets 
coherent states |5(, and squeezed vacuum pulses [1 
constitute an important step in demonstrating the poten- 
tial of this method. However, the efficiency and fidelity of 
the storage must be significantly improved before practi- 
cal applications become possible. 

In this paper, we present a comprehensive analy- 
sis of two recently demonstrated optimization proto- 
cols fl7l. 181 that are based on a recent theoretical pro- 
posallHllS EI US 13 ■ The first protocol iteratively 




optimizes the input pulse shape for any given control 
field while the second protocol uses optimal con- 
trol fields calculated for any given input pulse shape 
[HI]. We experimentally demonstrate their mutual con- 
sistency by showing that both protocols yield the same 
optimal control-signal pairs and memory efficiencies. We 
also show that for moderate optical depths (< 25), the 
experimental results presented here (as well as in Refs. 
[l7l HBl ) are in excellen t ag reement with a simple three- 
level theoretical model [l2l.Tl4j with no free parameters; 
we discuss the details of the correspondence between the 
actual atomic system and this simple model. Lastly, we 
study the dependence of memory efficiency on the optical 
depth. We show that for higher optical depths (> 25), 
the experimental efficiency falls below the theoretically 
predicted values; we discuss possible effects, such as spin- 
wave decay and four-wave mixing, that may limit the 
experimentally observed memory efficiency. 

The remainder of the paper is organized as follows: In 
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FIG. 1: (Color online) (a) The three-level A scheme used 
in theoretical calculations. The schematic (b) and example 
control (c) and signal (d) fields during light storage. At the 
writing stage (t < 0), an input signal pulse £i u (t) propagates 
through the atomic medium with low group velocity v g in the 
presence of a control field envelope Sl(t). While compressed 
inside the cell, the pulse is mapped onto a spin-wave S(z) by 
turning the control field off at time t — 0. After a storage 
period r, the spin-wave is mapped back into an output signal 
pulse £ ou t(t) using the retrieval control field envelope Q(t) 
(t > r). 



Sec. [TT1 we briefly summarize the three-level theory gov- 
erning the two procedures for optimizing photon storage. 
In Sec. IIII1 we describe our experimental system and dis- 
cuss its correspondence to the three-level model. In Sees. 
IIVI and [V] we present the results of experimental studies 
of both optimization procedures and demonstrate their 
consistency. In Sec. lVIl we investigate the dependence of 
memory efficiency on the optical depth of the medium. 
Finally, in Sec. IVII1 we conclude with the summary of 
our results. 
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II. REVIEW OF THE THEORY 

In this section, we briefly review the necessary con- 
cepts from the theoretical work [l2|, [H| on which our 
experiments rely. We consider the propagation of a weak 
signal pulse with envelope £(t) and a strong (classical) 
control field with a Rabi frequency envelope Q,(t) [19| 
in a resonant A-type atomic medium under the condi- 
tions of electromagnetically induced transparency (EIT), 
as shown in Fig. HJa). The control field creates a strong 
coupling between the signal field and a collective atomic 
spin excitation (spin wave) [4j. As a result, the initial 
pulse gets spatially compressed and slowed down inside 
the atomic ensemble. The group velocity of the pulse is 
proportional to the control field intensity (20j : 

Vg » 2\n\ 2 /( ai ) « c, (i) 

where 7 is the decay rate of the optical polarization and a 
is the absorption coefficient (i.e., unsaturated absorption 
per unit length), so that aL is the optical depth of an 
atomic medium of length L [2l[ . 

Fig. [ljb) illustrates schematically the three stages of 
the light storage process (writing, storage, and retrieval), 
while Figs. QJc) and[TJd) show control and signal fields, 
respectively, during a typical experimental run. At the 
writing stage, a signal pulse £ m (t) is mapped onto the 
collective spin excitation S(z) by adiabatically reducing 
the control field to zero. This spin wave is then preserved 
for some storage time r (storage stage), during which 
all optical fields are turned off. Finally, at the retrieval 
stage, the signal field £nut(t) is retrieved by turning the 
control field back on In the ideal case, the retrieved 

signal pulse is identical to the input pulse, provided the 
same constant control power is used at the writing and 
the retrieval stages. However, to realize this ideal storage, 
two conditions must be met. On the one hand, the group 
velocity v g of the signal pulse inside the medium has to 
be low enough to spatially compress the whole pulse into 
the length L of the ensemble and avoid "leaking" the 
front edge of the pulse past the atoms. This requires 
Tv g <C L, where T is the duration of the incoming signal 
pulse. On the other hand, all spectral components of 
the incoming pulse must fit inside the EIT transparency 
window to minimize spontaneous emission losses 1/T <C 
Aweit — V aLvg/L [3]. The simultaneous satisfaction 
of both conditions is possible only at very high optical 
depth aL > 1 

Experimental realization of very high optical depth in 
atomic ensembles requires high atomic density and/or 
large sample length. At high atomic density, EIT per- 
formance can be degraded by competing processes, such 
as stimulated Raman scattering and four-wave mixing 
[H H ill H HI S3 • Furthermore^spin-exchange col- 
lisions [28| | and radiation trapping [29|, [3(| HH may re- 
duce spin wave lifetime by orders of magnitude, limit- 
ing storage time and signal pulse durations. In addi- 
tion, achieving high optical depth in some experimental 



arrangements may be challenging, such as in magneto- 
optical traps (see e.g., Refs. @,|3|). Therefore, it is crucial 
to be able to maximize memory efficiency by balancing 
the absorptive and leakage losses at moderately large aL 
via optimal shaping of control and/or signal temporal 
profiles. To characterize our memory for light, we de- 
fine memory efficiency 77 as the probability of retrieving 
an incoming single photon after storage, or, equivalently, 
as the energy ratio between initial and retrieved signal 
pulses: 

/; +T |go Ut wi 2 ^ , 9 . 

V = n • \l) 

J° T \£ in (t)\ 2 dt 

The goal of any optimization procedure then is to max- 
imize r\ under the restrictions and limitations of a given 
system. 

In the theoretical treatment of the problem, the prop- 
agation of a signal pulse in an idealized three-level A sys- 
tem, shown in Fig. [IJa), is described by three complex, 
dependent variables, which are functions of time t and po- 
sition z [1, 0, • These variables are the slowly- varying 
envelope £ of the signal field, the optical polarization P 
of the \g) — |e) transition, and the spin coherence S. The 
equations of motion for these variables are [H, 0, [H| 

(d t +cd z )£(z,t) = igV~NP{z,t), (3) 
d t P{z,t) = -jP(z,t) + igVN£(z,t) + 

iCl(t- z/c)S(z,t), (4) 
d t S{z,t) = -j a S(z,t) + iCl(t- z/c)P(z,t){5) 

where g^/N = y^yac/2 is the coupling constant between 
the atomic ensemble and the signal field, and 7 and 7 S are 
the polarization decay rates for the transitions \g) — \e) 
and \g) — \s), respectively. While, in general, Eqs. 
cannot be fully solved analytically, they reveal several im- 
portant properties of the optimization process [13 . [l4T | . 
These properties are most evident in the case when spin 
wave decay rate 7 S is negligible during the processes of 
writing and retrieval (j s T < 1), which will hold for most 
of the discussion in the present paper, except for parts of 
Sec. IVII In this case, the highest achievable memory effi- 
ciency depends only on the optical depth aL and the mu- 
tual propagation direction of the control fields during the 
writing and retrieval stages [H[ . For each optical depth, 
there exists a unique spin wave, S op t{z), which provides 
the maximum memory efficiency. Thus, the focus of the 
optimization process becomes identifying a matching pair 
of writing control and signal pulses that maps the signal 
pulse onto this optimal spin wave. Note that no addi- 
tional optimization is required with respect to the re- 
trieval control field, because the memory efficiency does 
not depend on it, provided spin wave decay is negligible 
during retrieval [la [l4| . 

In the present experiments and in Refs. [13, [H| , the op- 
timization procedures are tested using weak classical sig- 
nal pulses rather than quantum fields. Such experimental 
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FIG. 2: (Color online) Experimental apparatus (see text for 
abbreviations). Inset: Schematic of the 87 Rb Di line level 
structure and relevant A systems formed by control and signal 
fields. 



arrangements greatly improved the experimental simplic- 
ity and the accuracy of data analysis. At the same time, 
the linear equations of motion for classical and quantum 
signal pulses are identical, which makes the presented re- 
sults applicable to quantized signal fields, such as, e.g., 
single photons. It is also important to note that the 
original theoretical work [l2], EH, EIL EE EH considered 
a wide range of interaction processes for storing and re- 
trieving photon wave packets (e.g., EIT, far-off-resonant 
Raman, and spin echo techniques) under a variety of con- 
ditions including ensembles enclosed in a cavity [13j . in- 
homogeneous broadening [l5l |. and high-bandwidth non- 
adiabatic storage (1/T ~ aLj) [li|. Since the proposed 
optimization procedures are, to a large degree, common 
to all interaction schemes and conditions, our results are 
relevant to a wide range of experimental systems. 



III. EXPERIMENTAL ARRANGEMENTS 

The schematic of the experimental apparatus is shown 
in Fig. [2] We used an external cavity diode laser (ECDL) 
tuned near the 87 Rb Di transition (A = 795 nm) with 
total available laser power w 45 mW. After separating a 
fraction of original light for a reference beam using a po- 
larizing beam splitter (PBS), the main laser beam passed 
through an electro-optical modulator (EOM), which 
modulated its phase at the frequency of the ground-state 
hyperfine splitting of 87 Rb (Ahf = 6.835 GHz) and pro- 
duced modulation sidebands separated by that frequency. 
We tuned the zeroth order (carrier frequency) field to the 
5 2 S 1/2 F = 2 -> 5 2 P 1/2 F' = 2 transition. This field was 
used as the control field during light storage. The +1 
modulation sideband played the role of the signal field 
and was tuned to the 5 2 S 1/2 F = 1 ->■ 5 2 P 1/2 F' = 2 
transition. 

To carry out the optimization procedure, we had to 



independently manipulate the amplitudes of the control 
and the signal fields. We used an acousto-optical modula- 
tor (AOM) to adjust the control field intensity. However, 
since all optical fields traversed the AOM, the intensities 
of all modulation comb fields were also changed. Thus, 
we accordingly adjusted the rf power at the EOM in- 
put (which controls the strength of the modulation side- 
bands) to compensate for any changes in the signal field 
amplitude caused by AOM modulation. 

To minimize the effects of resonant four-wave mix- 
ing, we filtered out the other (—1) first order modula- 
tion sideband (detuned by Ahf to the red from the car- 
rier frequency field) by reflecting the modulation comb 
off of a temperature-tunable Fabry-Perot etalon (FSR = 
20 GHz, finesse as 100). The etalon was tuned in reso- 
nance with this unwanted modulated sideband, so that 
most of this field was transmitted. At the same time, 
the control and signal field frequencies were far from the 
etalon resonance, and were reflected back with no losses. 
Such filtering allowed for suppression of the —1 modu- 
lated sideband intensity by a factor of « 10. 

Typical peak control field and signal field powers were 
18 mW and 50 fiW, respectively. The beam was weakly 
focused to w 5 mm diameter and circularly polarized 
with a quarter- wave plate (A/4). A cylindrical Pyrex 
cell (length and diameter were 75 mm and 22 mm, re- 
spectively) contained isotopically enriched 87 Rb and 30 
Torr Ne buffer gas, so that the pressure broadened optical 
transition linewidth was 27 = 2tt x 290 MHz [33] ■ The 
cell was mounted inside three-layer magnetic shielding 
to reduce stray magnetic fields. The temperature of the 
cell was controllably varied between 45° C and 75° C us- 
ing a bifilar resistive heater wound around the innermost 
magnetic shielding layer. 

We used relatively short pulses, so that spin decoher- 
ence had a negligible effect during writing and retrieval 
stages (except for parts of Sec. IVI[) and only caused a 
modest reduction of the efficiency during the storage time 
tx exp (— 27 s r). The Rb atom diffusion time out of the 
laser beam (~ 2 ms) was long enough to avoid diffusion- 
related effects on EIT dynamics [34 . |35| . We extracted 
the spin wave decoherence time by measuring the reduc- 
tion of the retrieved pulse energy as a function of storage 
time and fitting it to an exponential decay. We found the 
typical decay time to be 1/(2%) ~ 500 /is, most likely 
arising from small, uncompensated, remnant magnetic 
fields. 

After the cell, the output laser fields were recom- 
bined with the reference beam (at the unshifted laser 
frequency) at a fast photodetector, and the amplitude of 
each field was analyzed using a microwave spectrum ana- 
lyzer. Because of the 80 MHz frequency shift introduced 
by the AOM, the beatnote frequencies of the +1 and — 1 
modulation sidebands with the reference beam differed 
by 160 MHz, which allowed for independent measure- 
ment of the amplitude of each of these fields, as well as 
of the control field. 

To conclude this section, we explain the direct cor- 
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respondence between the experimental system and the 
theory based on three-level atoms [Fig. [lja)] that we re- 
viewed in Sec. HU The goal is to use the structure of 
the Di line of 87 Rb (see inset in Fig. [2]) to identify the 
optical depth aL and the control field Rabi frequency f2 
for the effective three level system. We first solve for the 
ground-state population distribution after control field 
optical pumping of the Rb D\ line, taking into account 
Doppler broadening, pressure broadening, and collisional 
depolarization of the excited state sublcvels 28]. We 
find the depolarization to be fast enough (for 30 Torr 
Ne, 7dopoi = 27r x 190 MHz (36|) to ensure roughly equal 
populations in each of \F,mp) = |1, — 1), |1, 0), |1, 1), 
and |2,2) ground state sublevels. Given this population 
distribution, we calculate the optical depth aL for the 
signal field as a function of Rb number density. For 
example, we find that at 60.5°C (Rb vapor density of 
2.5 x 10 11 cm" 3 ) the optical depth is aL = 24.0. More- 
over, approximately 60% of this optical depth comes 
from atomic population of \F,mp) = |1,1) due to the 
large corresponding Clebsch-Gordan coefficient. Thus, 
to calculate f2, we use the dipole matrix element of the 
\F = 2,m F = 1) -> \F' = 2,m F = 2) transition. 
Approximating a transverse Gaussian laser beam pro- 
file with a uniform cylindrical beam of diameter 5 mm 
of the same power, we find, for example, that for the 
control power of 16 mW, f2 = 2tt x 6.13 MHz. Since 
the collisionally broadened optical transition lincwidth 
(27 = 2n x 290 MHz) is comparable to the width of 
the Doppler profile, the effects of Doppler broadening 
are negligible, making Eqs. (0T5]) directly applicable. We 
note that all the theoretical modeling is done with no free 
parameters. 



IV. SIGNAL PULSE OPTIMIZATION 

One approach to the optimization of light storage is 
based on important time-reversal properties of photon 
storage that hold even in the presence of irreversible po- 
larization decay [Til ]. In particular, for co-propagating 
[32| writing and retrieval control fields, the following is 
true under optimized conditions (see Fig. [1]) : if a signal 
pulse £i n (t) is mapped onto a spin wave using a partic- 
ular control field f2(f) and retrieved after some storage 
time r using the time-reversed control field f2(T — t), 
the retrieved signal pulse shape £ ou t (t) is proportional to 
the time- reversed input signal pulse £i n (T — t), but at- 
tenuated due to imperfect memory efficiency. (Here and 
throughout the paper, control and signal envelopes are 
assumed to be real.) This symmetry also gives rise to 
an experimentally realizable iteration procedure, which, 
for any given writing control field, determines the op- 
timal incoming signal pulse shape. This procedure has 
been first demonstrated experimentally in Ref . [Tj] . The 
present experiment was performed independently on a 
different (although similar) experimental setup. There- 
fore, in order to use this procedure in Sec. lVII to study the 



£ OA- 
1 0.2- 
in 0.0 ■ 



-30 
Time ([is) 



■S OA- 

1 0.2- 
GO 0.0 ■ 

£0A- 
10.2 
000.0 ■ 

£0.4 

1 0.2 
in 0.0 



(1) 



-30 
Time (|js) 



effa 



-30 
Time ([is) 



- J 3 ' 4 ' 5 ) 
Bi„(t) 



£0.4- 
|0.2- 

00 0.0- 

£0A- 
1 0.2 
in 0.0 ■ 

£ OA- 

1 0.2- 
in 0.0 ■ 

£ OA- 
1 0.2- 
in 0.0 ■ 



Random initial pulse 



our ' 



'' I 1 ' 1 ' f i* ' 

50 100 
Time ([is) 



J- 



1 iteration 



-1— r—p- 1— 1— 1— 

50 1O0 
Time (us) 



J 



'■ ■'■ i 1 1 
50 ido 

Time (us) 



3 rd , 4 th and 5 ,h iterations 



-30 
Time (us) 



En,,,* 1 ) 



50 irjo 

Time (us) 



FIG. 3: (Color online) Iterative signal pulse optimization. 
The experimental data (solid black lines) is taken at 60.5 °C 
(aL — 24) using 16 mW constant control field during writing 
and retrieval (solid red line in the top panel) with a r = 100 /is 
storage interval. Numerical simulations are shown with blue 
dashed lines. Left: Input pulses for each iteration. Right: 
Signal field after the cell, showing leakage of the initial pulse 
for t < and the retrieved signal field £ out for t > 100 [is. 
Here and throughout the paper, all pulses are shown in the 
same scale, and all input pulses are normalized to have the 
J_ T \£i n (t)\ 2 dt = 1, where t is time in /is. 
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dependence of memory efficiency on the optical depth, 
we verify in this section its successful performance in the 
present experimental setup. In addition, the implementa- 
tion of iterative signal optimization in this experimental 
setup will allow us, in Sec. [Vj to compare and verify the 
consistency of signal and control optimizations. 

The sequence of experimental steps for the iterative 
optimization procedure is shown in Fig. [3] The plots 
show the control field and the measured and simulated 
signal fields (solid red lines in the top panel, solid black 
lines, and dashed blue lines, respectively). Before each 
iteration, we optically pumped all atoms into the state 
\g) by applying a strong control field. We started the op- 
timization sequence by sending an arbitrary signal pulse 
£jj/(t) into the cell and storing it using a chosen con- 
trol field f2(t). In the particular case shown in Fig. O 
the group velocity was too high, and most of the in- 
put pulse escaped the cell before the control field was 
reduced to zero. However, a fraction of the pulse, cap- 
tured in the form of a spin wave, was stored for a time 
period r = 100 /is. We then retrieved the excitation 
using a time-reversed control field fl(t) = f2(r — t) and 

recorded the output pulse shape £^ t (t). For the sam- 
ple sequence shown, the control fields at the writing and 
retrieval stages were constant and identical. This com- 
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FIG. 4: (Color online) (a) Experimental (solid) and theo- 
retical (dashed) optimized signal pulses obtained after five 
steps of the iteration procedure for three different powers of 
the constant control fields during writing and retrieval stages, 
(b) Corresponding memory efficiencies determined for each it- 
eration step. Theoretically predicted optimal efficiency value 
is shown by the dashed line. The temperature of the cell was 
60.5 °C (aL = 24). 



pletes the initial (zeroth) iteration step. The efficiency 
of light storage at this step was generally low, and the 
shape of the output pulse was quite different from the 
time-reverse of the initial pulse. To create the input 
pulse (t) for the next iteration step, we digitally time- 
reversed the output foutW °f the zeroth iteration and 
renormalized it to compensate for energy losses during 
the zeroth iteration: £^' (t) oc £^ t (r — t). Then, these 
steps were repeated iteratively until the rescaled output 
signal pulse became identical to the time-reversed profile 
of the input pulse. As expected, the memory efficiency 
grew with each iteration and converged to 43 ± 2%. 

To verify that the obtained efficiency is indeed the 
maximum possible at this optical depth and to confirm 
the validity of our interpretation of the results, we com- 
pare the experimental data to numerical simulations in 
Fig. O Using the calculated optical depth and the con- 
trol Rabi frequency (see Sec. IIII[) , we solve Eqs. ([3][5]) 
analytically in the adiabatic limit TaL^f 1 14|, which 
holds throughout this paper. There is a clear agreement 
between the calculated and measured lineshapes and am- 
plitudes of the signal pulses. Also, theory and experiment 
converge to the optimal signal pulse shape in a similar 
number of iteration steps (2-3), and the experimental ef- 
ficiency (43 ± 2%) converged to a value close to the the- 
oretical limit of 45% (see below). 

As in our previous study [l7| . we confirmed that the 
final memory efficiency and the final signal pulse after a 
few iteration steps are independent of the initial signal 

pulse £^\t). We also confirmed that the optimization 
procedure yields the same memory efficiency for differ- 
ent control fields. While constant control fields of three 
different powers yield different optimal signal pulses [Fig. 
SJa)], the measured efficiency [Fig. Hfb)] converged after 
a few iteration steps to the same value of 43 ± 2%. With 
no spin wave decay, the highest achievable memory effi- 
ciency for the optical depth aL — 24 is 54% [l4| ■ Taking 
into account spin wave decay during the 100 fj,s storage 



time by a factor of exp[— 100/zs/500/is] = 0.82, the high- 
est expected efficiency is 45% [dashed line in Fig. H£b)], 
which matches our experimental results reasonably well. 



V. CONTROL PULSE OPTIMIZATION 

The iterative optimization procedure described in the 
previous section has an obvious advantage: the optimal 
signal pulse shape is found directly through experimental 
measurements without any prior knowledge of the sys- 
tem parameters {e.g., optical depth, control field Rabi 
frequency, various decoherence rates, etc.). However, in 
some situations, it is difficult or impossible to shape the 
input signal pulse (e.g., if it is generated by parametric 
down-conversion [37j|). In these cases, the control field 
temporal profile must be adjusted in order to optimally 
store and retrieve a given signal pulse. 

To find the optimal writing control field for a given in- 
put pulse shape £{ n (t), we maximize r\ [Eq. ^] within the 
three-level model [Eqs. (JHSSJ)] . In this model, for a given 
optical depth aL and a given retrieval direction (coincid- 
ing with the storage direction in the present experiment 
(32j|). there exists an optimal spin wave S op t(z), which 
gives the maximum memory efficiency. One way to cal- 
culate the control field required to map the input pulse 
onto this optimal spin wave is to first calculate an artifi- 
cial "decayless" spin wave mode s(z), which, like S op t(z), 
depends only on the optical depth and not on the shape 
of the incoming pulse. This "decayless" mode s(z) hy- 
pothetically allows for unitary reversible storage of an 
arbitrary signal pulse in a semi-infinite and polarization- 
decay-free atomic ensemble, in which the group velocity 
of the pulse is still given by Eq. ([1]). The unitarity of the 
mapping establishes a 1-to-l correspondence between a 
given input signal pulse shape £i n (t) and an optimal writ- 
ing control field that maps this input pulse onto s(z). The 
same control field maps this input pulse onto the true op- 
timal spin wave S op t(z), once polarization decay and the 
finite length of the medium are taken into account. The 
details of this construction are described in Ref. [l4l |. 

As an example of control field optimization, we con- 
sider the storage of three different initial pulse shapes, 
shown by dotted black lines in the middle row in Fig. [5] 
a step with a rounded leading edge (a'), a segment of 
the sinc-function (b'), and a descending ramp (c'). The 
top row (a,b,c) shows the corresponding calculated op- 
timal writing (t < 0) control pulses. Since the shape 
and power of the retrieval control pulse do not affect the 
memory efficiency fl2l fl4| . we show, in the top row of 
Fig. [5l two retrieval control fields for each input pulse: 
a flat control field (dashed) and the time-reverse of the 
writing control (solid). As expected, the flat control field 
(the same for all three inputs) results in the same out- 
put pulse [dashed in (a', b', c')] independent of the input 
signal pulse, because the excitation is stored in the same 
optimal spin wave in each case. On the other hand, us- 
ing the time-reversed writing control field for retrieval 
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FIG. 5: (Color online) Storage of three signal pulses (a', b', c') using calculated optimal storage (t < 0) control fields (a), (b), 
(c). Input signal pulse shapes are shown in black dotted lines. The same graphs also show the leakage of the pulses (solid black 
lines for t < 0) and retrieved signal pulses (t > 100 /is) using flat control fields at the retrieval stage (dashed red lines), or using 
time-reversed control fields (solid red lines). Graphs (a", b", c") show the results of numerical calculations of (a', b', c'). The 
temperature of the cell was 60.5 °C (aL = 24). 



yields output pulses that are time-reversed (and atten- 
uated) copies of the corresponding input pulses. This 
means that the time-reversal iterations of Sec. IIVI start- 
ing with these control-signal pairs converge on the zeroth 
iteration, which proves the consistency of the signal opti- 
mization of Sec. IIVI with the control optimization of the 
present section. The experimental data also agrees very 
well with numerical simulations [bottom row (a", b", c") 
in Fig. [5], supporting the validity of our interpretation of 
the data. 

To further test the effectiveness of the control opti- 
mization procedure, we repeated the same measurements 
for eight different randomly selected pulse shapes, shown 
as black lines in Fig. [HJa). Pulses #4, #6, and #8 are 
the same as the input pulses (a'), (b'), and (c') in Fig. [5] 
For each of the eight input pulses, we calculated the opti- 
mal writing control [red lines in Fig.^a)] and then mea- 
sured the memory efficiency [Fig. [SJb)], retrieving with 
either a constant control pulse or a time-reversed writing 
control pulse (open red diamonds and solid black cir- 
cles, respectively). The measured efficiencies are in good 
agreement with each other and with the theoretically cal- 
culated maximum achievable memory efficiency of 45% 
(horizontal dashed line) for the given optical depth. 

By performing these experiments, we found that 
knowledge of accurate values for the experimental param- 
eters, such as optical depth or control field intensity, is 
critical for calculations of the optimal control field. Even 
a few percent deviation in their values caused measur- 
able decreases in the output signal pulse amplitude. In 
our experiment, effective optical depth and control field 
Rabi frequency were computed accurately directly from 
measurable experimental quantities with no free parame- 



ters. The accuracy of the parameters was also verified by 
the excellent agreement of experimental and theoretical 
results of iterative optimization in Sec. IIVI We note that 
for some other systems, the necessary experimental pa- 
rameters may be difficult to compute directly with high 
accuracy; in that case, they can be extracted from the 
iteration procedure of Sec. |TVl 



VI. DEPENDENCE OF MEMORY EFFICIENCY 
ON THE OPTICAL DEPTH 

In the previous two sections, we verified at optical 
depth aL = 24, the consistency of the signal and con- 
trol optimization methods and their agreement with the 
three-level theory. In this section, we study the depen- 
dence of memory efficiency on optical depth. To verify 
the theoretical prediction that the optimal efficiency de- 
pends only on the optical depth of the sample, we re- 
peated the iterative signal optimization procedure (Sec. 
IIVP for several constant control field powers at different 
temperatures of the Rb cell ranging from 45°C (aL = 6) 
to 77°C (aL = 88). In Fig.[7fa), we plot the measured ef- 
ficiencies (markers) along with the maximum achievable 
efficiency predicated by the theory without spin decay 
(thin black line) and with spin decay during the storage 
time (thick black line). This graph allows us to make 
several important conclusions. 

First of all, it demonstrates that for relatively low op- 
tical depths (aL < 25), the optimized memory efficiency 
for different control fields is the same, to within the ex- 
perimental uncertainty, and approximately matches the 
theoretical value (thick black line). This confirms that 
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FIG. 6: (Color online) (a) Eight randomly selected signal pulse shapes (black lines) and their corresponding optimal control 
fields (red lines), (b) Memory efficiency for the eight signal pulse shapes using calculated optimized control fields at the 
writing stage, and flat control fields (open red diamonds) or inverted writing control fields (solid black circles) at the retrieval 
stage. Theoretically predicted optimal memory efficiency is shown by a dashed line. The temperature of the cell was 60.5 °C 
(aL = 24). 
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FIG. 7: (Color online) Memory efficiency as a function of optical depth obtained by carrying out iterative signal optimization 
until convergence, (a) At each optical depth, we considered constant control fields at four different power levels (indicated on 
the graph) during writing and retrieval stages. Note that many experimental data points overlap since the converged efficiencies 
are often the same for different control fields. Dashed lines are to guide the eye. Thin and thick black solid lines show the 
theoretically predicted maximum efficiency assuming no spin-wave decay and assuming an efficiency reduction by a factor of 
0.82 during the 100 /j,s storage period, respectively, (b) Thin and thick black lines are the same as in (a), while the three lines 
with markers are calculated efficiencies for three different control fields (indicated on the graph) assuming spin wave decay with 
a 500 us time constant during all three stages of the storage process (writing, storage, retrieval). 



the optimization procedure yields the maximum effi- 
ciency achievable for a given optical depth. However, 
for aL > 20, the efficiency obtained with the lowest con- 
trol field power (black empty circles) dropped below the 
efficiency obtained for higher control powers. As we will 
now show, the most probable reason for such deviation 
is spin wave decay during writing and retrieval. 

As the optical depth increases, the duration of the op- 
timal input pulse increases as well, as shown in Fig. [SJa), 
following the underlying decrease of group velocity: T ~ 
L/vg oc aL Thus, above a certain value of aL, the 
duration of the optimal pulse for a given control field 
becomes comparable with the spin wave lifetime, and 
the spin wave decoherence during storage and retrieval 
stages can no longer be ignored. Further increase of the 
optical depth leads to a reduction of retrieval efficiency, 
even though the iterative optimization procedure is still 



valid [14] and produces signal pulses that are stored and 
retrieved with the highest efficiency possible for a given 
control field and aL. Fig.[7][b) shows the calculated max- 
imum achievable efficiencies for different constant control 
powers as a function of the optical depth, taking into ac- 
count spin wave decay with a 500 /zs time constant dur- 
ing all three stages of light storage. For each control field 
power, the efficiency peaks at a certain optical depth, and 
then starts to decrease as optical depth increases further. 
Since lower control powers require longer optimal input 
pulses T ~ L/v g oc 1/|^| 2 [see Fig.|U[a)], the correspond- 
ing efficiency reaches its peak at lower optical depths. 
Thus, the problem of efficiency reduction posed by spin- 
wave decay during writing and retrieval can be alleviated 
by using higher control powers, and hence shorter opti- 
mal signal pulses. While this effect explains the reduction 
of maximum memory efficiency attained with the lowest 
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FIG. 8: (Color online) Results of the optimization procedures for different optical depths: aL = 24 (red), aL — 40 (black), and 
aL — 50 (green). The top panel [(a) and (b)] shows storage and retrieval (b) of the optimized input signal pulses (a) obtained 
by running iterative optimization until convergence for a constant control field of power 8mW [dash-dotted line in (b)] . Solid 
lines correspond to experimental results, while dashed lines show the results of numerical simulations. In the bottom panel 
[(c) and (d)], (c) shows the calculated optimal writing control fields (t < 0) for a step- like signal pulse [dotted line in (d)] and 
the time-reverses of these control fields used during retrieval (t > 100 /is), while (d) shows the resulting storage followed by 
retrieval. 



control power for aL > 20 [Fig. [2a)], other effects, dis- 
cussed below, degrade the efficiency for all other control 
powers for aL > 25, as indicated by the divergence of 
experimental data in Fig. [7Ja) from the corresponding 
theoretical efficiencies in Fig. [Tfb) (red and green lines). 
Remarkably, at these optical depths, the iterative signal 
optimization procedure still yields efficiencies that grow 
monotonically at each iteration step for the three highest 
control powers. This suggests that iterative signal opti- 
mization may still be yielding the optimum efficiency, 
although this optimum is lower than what the simple 
theoretical model predicts. 

To further test the applicability of our optimization 
procedures at higher optical depths, we complemented 
the signal-pulse optimization [Fig. [H]Ja,b)] with the cor- 
responding control field optimization [Fig. [8jc,d)]. We 
stored and retrieved input pulse #4 from Fig. [Hfa) us- 
ing calculated optimal writing control fields [t < in 
Fig. [5)Jc)] at different optical depths aL — 24, 40, and 
50. As expected, the overall control power was higher 
at higher optical depths to keep the group velocity un- 
changed: L/T ~ v g (x Vi 1 /(aL). For each optical depth, 
we used a time-reversed writing control field to retrieve 
the stored spin wave. This resulted in the output signal 
pulse shape identical to the time-reversed (and atten- 
uated) copy of the input pulse, as shown in Fig. [8jd). 
Although the memory efficiency drops below the theo- 
retical value at these high optical depths [aL = 50 for 
the green lines in Fig. |8Jc,d)], the results suggest that 
the calculated control field may still be optimal, since it 
yields the time- reverse of the input signal at the output. 



To gain insight into what may limit the memory ef- 
ficiency for 25 < aL < 60, we investigated the effect of 
resonant four-wave mixing. Thus far, we have considered 
only the ground-state coherence created by the control 
and signal fields in the one-photon resonant A configura- 
tion [Fig.[]Ja)]. However, the strong control field applied 
to the ground state \g) can also generate an additional 
Stokes field Eg, as shown in Fig.^a). This process is sig- 
nificantly enhanced in EIT media [22|, . In particular, 
it has been shown that a weak signal pulse traversing an 
atomic ensemble with reduced group velocity generates a 
complimentary Stokes pulse that travels alongside with 
a comparably low group velocity (38l . l39j . 

To determine the effect of resonant four-wave mixing 
on light storage, we first carried out iterative signal op- 
timization for a constant control field pulse of 16 mW 
power at different optical depths, but then detected not 
only the signal field, but also the Stokes field, at the 
retrieval stage [see Fig. [3£b)]. We see that at low opti- 
cal depths, the retrieved Stokes pulse [blue empty dia- 
monds] is negligible compared to the output signal pulse 
[red filled diamonds, which are the same as the red filled 
diamonds in Fig. EJa)]. However, at aL > 25, the energy 
of the output pulse in the Stokes channel becomes sig- 
nificant. While the energy of the retrieved signal pulse 
stayed roughly unchanged for 25 < aL < 60, the energy 
of the output Stokes pulse showed steady growth with 
increasing aL. Moreover, the combined energy (black 
empty circles) of the two pulses retrieved in the signal 
and Stokes channels added up to match well the theoret- 
ically predicted highest achievable efficiency (solid black 
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FIG. 9: (Color online) (a) Level diagram illustrating Stokes 
field (Es) generation due to resonant four- wave mixing, (b) 
Memory efficiency for retrieval in the signal channel [same as 
the red filled diamonds in Fig. E£a)], Stokes channel, and the 
total for both channels. The efficiencies are obtained by car- 
rying out iterative optimization till convergence for constant 
writing and retrieval control fields of 16 mW power. Dashed 
lines are to guide the eye. The solid line (same as the thick 
black line in Fig. [7]) shows the theoretically predicted maxi- 
mum efficiency assuming an efficiency reduction by a factor 
of 0.82 during the 100 /is storage period. 



line). We will study elsewhere whether this match is 
incidental and whether it can be harnessed for memory 
applications. For the purposes of the present work, we 
simply conclude that the effects of four-wave mixing be- 
come significant around the same value of aL (~ 25) 
where experiment starts deviating from theory. There- 
fore, four-wave mixing may be one of the factors respon- 
sible for the low experimental efficiencies at high optical 
depths. 

For aL > 60, iterative signal optimization still con- 
verges, but efficiency does not grow monotonically at 
each iteration step, which clearly indicates the break- 
down of time-reversal-based optimization [l4l ]. In addi- 
tion, the final efficiency is significantly lower than the 
theoretical value (Fig.[7|). Many factors, other than four- 
wave mixing, may be contributing to the breakdown of 



time-reversal-based optimization and to the rapid de- 
crease of memory efficiency at aL > 60. First of all, the 
absorption of the control field at such high optical depths 
is significant (measured to be > 50%). In that case, the 
reabsorption of spontaneous radiation contributes appre- 
ciably to spin wave decoherence [29l . [30l | and can make 
the spin wave decay rate 7 S grow with aL, reducing the 
light storage efficiency [3l| . Spin-exchange collision rate 
[28j |. which destroys the spin- wave coherence, also be- 
comes significant at high Rb density, reducing spin wave 
lifetime even further. 



VII. CONCLUSIONS 

We have studied in detail two quantum memory opti- 
mization protocols in warm Rb vapor and demonstrated 
their consistency for maximizing memory efficiency. We 
have also observed good agreement between our experi- 
mental data and theoretical predictions for relatively low 
optical depths (< 25), both in terms of the highest mem- 
ory efficiency and in terms of the optimized pulse shapes. 
At higher optical depths, however, the experimental ef- 
ficiency was lower than predicted. We observed that 
resonant four-wave mixing processes became important 
at these higher optical depths. We expect our studies 
to be of importance for enhancing the performance of 
ensemble-based quantum memories for light. 
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